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Abstract — We propose theoretically a generalized memristive 
system based on controlled spin polarizations in the giant 
magnetoresistive material using a feedback loop with the classical 
Hall Effect. The dynamics can exhibit a memristive pinched 
hysteretic loop that possesses the self-crossing knot not located 
at the origin. Additionally, one can also observe a single-looped 
orbit in the device. We also provide a sufficient condition for 
the stability based on an estimation of the Floquet exponent. 
The analysis shows that the non-origin-crossing dynamics is 
generally permitted in a class of passive memory systems that 
are not subject to Ohm's Law. We further develop the prevailing 
homogeneous definition to a broadened concept of generalized 
heterogeneous memristive systems, permitting no self-crossing 
knot at the origin, and ultimately to the compound memory 
electronic systems. 

Index Terms — memristor, nonlinear dynamics, spintronics, 
Hall Effect. 

I. Introduction 

IN 1971, the concept of the memristor was coined by Leon 
O. Chua at Berkeley Q. But it was not until 2008 that 
S. Williams and his group at Hewlett-Packard Laboratories 
shed the first light on this scientific preemie with their 
T1O2-X memristor Q. It was also from then on that nonlinear 
electronics has been recognized as one of the comprising 
bases of the entire mansion of electronics. No later than this 
revival had Di Ventra, et al, pushed the concept forward 
to a large family of passive memory systems, consisting of 
memristive, memcapacitive and meminductive systems 0. 
Deeper insights have also permeated into the territory of 
biophysiology, building links between neuronal transmissions, 
learning rules and memristive behaviours (4)- (7), making it 
an indeed interdisciplinary research subject. Parallel to the 
conceptual development, researchers have been enthused in 
seeking memristive systems on new physical mechanisms, at 
both the material level and the device level ifHTl - iTTsTl . Among 
them is the spintronics of great interest, which has linked 
the memristive concept to amounts of established research 
results in this area. The spin Hall Effect in a semicon- 
ducting system with an inhomogeneous electron density has 
demonstrated a hysteretic pinch under periodically driving 
voltages fT2ll . and the spin diffusion and relaxation processes 
in the semiconductor/half-metal junction can also exhibit a 
memristive behaviour fT3ll . On the other hand, in the metal 
material, memristive systems are realizable by employing spin- 
torque-induced magnetization switching or magnetic-domain- 
wall motion (HIGH- While people have been focused on 
other spin polarization mechanisms, in this paper we propose 
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theoretically a memristive system based on the controlled 
spin polarization in the giant magnetoresistive material using 
a feedback loop with the classical Hall Effect. The system 
is of special interest in exhibiting hysteretic-pinched (as in 
a standard memristor) loop, which, however, has the self- 
crossing knot not at the origin. This characteristic has grown 
out of the prevailing categorization of memristive systems 
proposed by Di Ventra, et al, but is still within a more general 
frame of memristive systems, as also suggested by Chua 
recently. It is hence defined as a generalized memristor model, 
which we name as the heterogeneous memristive system. 
The paper is organized as follows: we will introduce the 
generalized memristive systeml in Sec. II. The memristive 
dynamics is depicted and analyzed in Sec. Ill, including a 
sufficient condition for stability. In Sec. IV, we will relocate at 
the system's quasi- standard memristive dynamics, and discuss 
the generalization to the concept of heterogeneous memristive 
systems and compound memory systems. We will draw the 
conclusions in Sec. V. 

II. The Proposed Spin-Valve Memristive System 

We propose a spin-valve memristive system established on 
a feedback loop. The loop controls the spin polarization in a 
giant magnetoresistive material via the classical Hall Effect, as 
shown in Fig. [T] which gives rise to a memory behaviour. Let 
us first consider a giant magnetoresistive (GMR) material in 
the form of stacked thin-films based on the spin- valve mech- 
anism, of which the resistance reveals a strong dependence 
on the applied magnetic field lfT6ll[T8l . The principle of the 
significant change in the resistance is caused by the parallel 
or antiparallel spin polarizations: an electron passing through 
the GMR material will be scattered more if the spin of the 
electron is opposite to the direction of the magnetisation in the 
ferromagnetic layer, and the material hence expresses a larger 
resistance r^; otherwise, it expresses a smaller resistance r^. 
With this mechanism, an equivalent circuit for a stack of films 
can be constructed, so that the compound resistance under a 
magnetic field (corresponding to parallel spin alignment 
in the material) is much lower than that without the magnetic 
field R^i (corresponding to unparallel spin alignment in the 
material). For our interest in the study of the dynamics, we 
can characterize the resistance empirically by a hyperbolic 
function of the magnetic flux <p m as 

R(^m) = Rn + («u - R n) sech (f^) (D 

\Pm0j 

with (/> m0 denoting a normalization quantity with the dimen- 
sion of magnetic flux. The hyperbolic function is a fitting 
curve to the characteristic curves of a typical class of GMR 
materials. For example, as shown in Fig. [T] the hyperbolic 
curve (dashed line) given in Eq.[T]can well fit the experimental 
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results (solid line) in the report of A. Fert, et.al, where the 
thin-film structure is composed of alternating ferromagnetic 
and non-magnetic layers, namely, ... |Fe|Cr|Fe|Cr| with the 
thickness of Fe at 3 nm and Cr at 0.9^1.8 nm lfT6l . On 
the other hand, the sech function has a parabolic form when 
expanded to the second order for small fields, which has a 
decent match with the existing theorectical formalizations (for 
independent moment models) in |[T9l and l20l . In spite of 
various GMR structures with specific quantum models, this 
empirical formalization deals with the problem of interest 
without losing generality. The characteristic curves of GMR 
materials are also not limited to the form given in fT6lL but 
they will generally give similar memristive dynamics in our 
proposed structure. 

To give rise to a memristive effect, we employ the classical 
Hall Effect to feed back a magnetic field as a control of 
the spin polarization in the ferromagnetic layers of the GMR 
material (the red-coloured layers in Fig. [T]). The classical Hall 
Effect is based on a balance between the Lorenz force and the 
established inner electric field l2T1l . The Hall electromotive 
force Eh is related to a magnetic field on the conductor by 
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where n is the carrier density in the magnetoresistive conduc- 
tor, e is the absolute value of the electron charge, d is the 
transverse scale of the conductor, S = ab is the horizontal 
cross-section area of the conductor. However, in this model, 
the magnetic field B is supplied internally by the Hall current 
in- When the Hall current is drawn by the Hall electromotive 
force out to a coil, feeding back a controlling magnetic flux m 
on the conductor itself, the balance is sustained in a dynamic 
way. The electromotive force equals the voltage drop on the 
coil 



£h = 



d^rn 

dt 
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where Rl is the coil resistance. It is to note that attributed to 
the ferromagnetic core, the inductance L increases to a large 
extent with the enhanced permeability (by a factor of several 
hundred to thousand) and the dimension of the coil is allowed 
to be larger than that of the core due to the magnetic flux 
confinement. On the other hand, with the energy conservation, 
the entire system absorbs the power offered by an external 
voltage u(t), equaling to the dissipation in the resistance of 
the conductor and the power in the coil: 
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in which it is especially notable that due to the extraction 
of the Hall current, Ohm's law does not simply apply. In 
order to reveal the dominant factors contributing to the mem- 
ristive dynamics, we base our modelling on the following 
approximations: l.the magnetic flux changes with the Hall 
current in linearly through an inductance L of the coil, i. 

(j) m = L • in. However, in terms of accuracy, it is worth 
pointing out that the ferromagnetic core saturation can cause a 
deviation from the linear relation at a high coil current. 2. the 
coil resistance Rl is constant. In a more accurate manner, 
this resistance also includes a transverse resistance of the 
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Fig. 1. The proposed feedback spin-valve memristor model in a giant magne- 
toresistive material of the current-perpendicular-to-plane (CPP) configuration, 
(a) structure: the spin in the ferromagnetic layers (red) are polarized adaptively 
to the temporally changing magnetic flux induced by the Hall current, (b) the 
principle of spin- valve in the GMR material and the equivalent circuit model, 
(c) the relation between the giant magnetoresistance R(4>m) and the magnetic 
flux 4>m described by an empirical formula with a hyperbolic function, which 
is here used as a fitting curve to the experimental result in 1161 . The solid 
line is the experimental curve, and the dashed line is the fitting hyperbolic 
curve. 
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ferromagnetic conductor. For simplicity, we suppose that the 
transverse resistance is relatively small (by properly choosing 
the size ratio of the GMR material) and its temporal change 
is negligible in series with the main coil resistance, and hence 
Rl is regarded constant. 3. other energy losses in the GMR 
material are not counted in for simplicity, which are mainly 
caused by the Eddy currents and the coercive losses (see, e. 
g., the GMR material in Ref. |[T6l with a low coercitivity). 
These approximations aid to simplify the model, while also 
maintains the factors determining the dynamics, esp. the 
orbital topologies of our main interest. The nonlinearity in the 
inductance of the coil will limit the magnetic flux in practice. 
When these non-ideal effects are involved, they will deform the 
orbital shapes and modify the parameters into a more realistic 
range, but the system will not lose the memristive dynamics in 
general. Based on the above factors and approximations, the 
complete dynamics of the proposed memristive system can 
eventually be defined by the following equation set: 

U(t) R{<t> m )i + ^-(t> 2 m 



(5) 



with 



and 



R^m) = Rft + (R 
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dt -^" L " L 
in which Dq denotes the coefficient 1/nedS. It is supposed 
here that the average exertion time of spin torques under 
the magnetic field is sufficiently short compared to the time 
scale of the dynamics, and hence no time delay is involved 
in the above model. Special notice is to take on the extra 
term deviating to the resistive relation in Eq. ([5]), which is 
responsible for giving rise to a new memory effect beyond 
standard memristive systems. For material implementation, 
the current-perpendicular-to-plane (CPP) configuration of thin- 
films is to adopt (as shown in Fig. [TJ, for yielding the highest 
giant magnetoresistance along the external current l22l - l24l . 
which possesses a coefficient GMR = (R^ — R^)/R^ 
exceeding 1 at the present time. 

III. The Dynamics of the Memristive System 

The numerical calculations display that the proposed system 
can exhibit a memristive hysteretic pinch in the voltage- 
current phase plane when driven by a sinusoidal voltage signal 
(u(t) = uo cos(uit)), as shown in Fig. Qa)), which realizes 
a memory effect of the dynamics produced by the feedback 
loop. However, unlike the typical hysteretic pinch in standard 
memristors, this generalized memristive system has its self- 
crossing knot of the orbit deviated from the origin of the 
phase plane, though the deviation can be rather small. This 
interesting characteristic is caused by the heterogeneous term 
in Eq. ([5]), in which the voltage u is no loner a homogeneous 
equation about the current i with the L-containing term. In 
this regarding, we will name such a system as a Heteroge- 
neous Memristive System, as a generalization of the standard 
memristive system, which we will name specifically as the 
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Fig. 2. The memristive dynamics of the spin- valve system, with (a) the 
sinusoidal input signal and (b) the square-wave input signal (hard switching 
response). The parameters are chosen as e = 1.6 x 10 -19 C, n = 
10 20 /cm 3 ,i? n = 0.33 Q,R n = 0.66 Q, Rl = 0.5 ft, a = 1 fim,b = 
6 //in, d = 1 /yjn, L = 50 nH and (f) m o = 30 pWb. The sinusoidal 
and square- wave amplitudes are both uo = 80 mV at the frequency 
/ = uo /27T = 10 GHz. The zoomed-in display the deviation of the trajectory 
from the origin. The initial condition is set to be </> m (0) = 9 pWb. 



Homogeneous Memristive System. We will discuss in details 
the significance of this categorization in Sec. IV. For this 
specific system, it is worth pointing out that the proposed 
memristive system can be regarded as a standard memristive 
system in a limit. In fact, if we simplify Eq. ([5]) and ([6]) 
by diminishing the terms containing L when these terms are 
small- valued (e.g., with a relatively large L) and integrating 
Eq. ([6]), we can recognize this degenerate model as a charge- 
controlled memristor with the memristance M(q) = R((j> m ), 
by taking </>^ as the state- variable: 



u{t) = R(y/fa)i 
dln((f) 2 m /(f) 2 m0 ) 



dt 



= 2D i 



(7a) 
(7b) 
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noting that R((j)m) is an even function of m . Henceforth, 
the proposed memristive system can be regarded as charge- 
controlled memristor in this limit. Hard switching also demon- 
strates a hysteretic loop, as shown in Fig. [2jb). 

Regarding the local properties, we use the instantaneous 
differential resistance Rd(t) = du/di = (du/dt) • (di/di) -1 
to characterize the orbits, as shown in Fig. [3] The behavior of 
Rd(t) is a compound expression of both the dynamical process 
of the giant magnetoresistance and the non-Ohmic effect. 
We divide the Rd(t) into two sections: the upward section 
R^(t) = Rd(t | u(t) : —uq —> +^o) and the downward 
section Rd(t) = Rd(t | u(t) : -f^o —> —v>o)- The differ- 
ential resistance of the upward section at the orgin R^(t ) 
approximates to (when u(t Q ) = and d(j> m {t ) / dt « 0, 
(j> m (t ) ~ 0), while the downward section R^(t) approaches 
R^ at the self-crossing knot t c . 

We have also observed that the system can also evolve to 
a memory system, which expresses a single-looped orbit (see 
Fig. Qb)). Such orbit rises up when the L-containing terms 
gain dominancy in Eq. §5§ and (b). In contrast to the standard 
memristive systems, the single-looped orbit also has a part of 
the trajectory that does not cross the origin when the voltage is 
zero. This is also simply attributed to the additional non-zero 
term to the resistive relation in Eq. ([5]). The hard switching 
also demonstrates a single-looped orbit, as shown in Fig.|5Jc). 

Globally, the bi-looped and single-looped orbits are topo- 
logical^ homeomorphic to the Lissajous figures for the ratio 
of 1 and 2, respectively. In fact, if we look at the dynamics in 
the u — <j) m phase plane, the two orbits are both homeomorphic 
to a circle. The uprising of the single-looped orbit is because 
the terms containing L gain dominancy in Eq. ^ and (b). 

We are interested now in the condition for stability. Given 
proper parametric conditions, this proposed memristive sys- 
tems can display a stable behavior after experiencing a tran- 
sient. It is observable that the dynamics of the model in the 
u — <pm phase plane displays a one-looped limit cycle under a 
sinusoidal driving signal, e. g., as displayed in Fig. [5] But 
to confirm it as an attractor, it is necessary to know the 
parameter range for stability. The model has one state-variable, 



Fig. 4. The single-looped dynamics of the spin- valve system, with a 
sinusoidal input signal. The parameters are chosen ase = 1.6xl0 -19 C, n = 
10 20 /cm 3 ,i? n = 0.33 fl,R n = 0.66 Q,R L = 0.5 ft, a = 1 fjm,b = 
6 /im, d = 1 /im, L = 4 nH and m o = 5 pWb. The sinusoidal amplitudes 
is uo = 20 mV at the frequency / = uj/2tt = 0.25 GHz. The zoomed-in 
display the deviation of the trajectory from the origin. The initial condition 
is set to be </> m (0) = 9 pWb. 



as rewritten in the following form 

d(j) r n D 



dt 

u{t) = uq cos(cjt) 



(j> m u{t) 
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We here use the Poincare map to study the stability of the 
periodic u — (j) m planar orbit [25]. The method converts 
a problem about a closed orbit into a problem about the 
fixed points of a mapping. However, the rare possibility of 
achieving an analytical expression of the mapping constrains 
most of on-hand analysis on numerical calculations for specific 
parameters. In order to see the parameter range for stability, 
we here rather provide a sufficient condition for the asymptotic 
stability by estimating the Floquet exponent of the system. 

Let 0^(t) be the closed orbit and rj(t) be an infinitesimal 
perturbation: m (t) = <j>m{t) +Tj(t). Substituting this relation 
into Eq. ( [8a| and using the Taylor expansion to the first order 
of 77, we can rewrite the equation in the following form: 



d((j)*m + rj) 
dt 



{ \d (R u - Rn) x^ sech ( J tanh 

L L 0mO \<fim0 J \(Pm0 
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Since 4>m(t) i s a solution of Eq. (8a), the following relation 
should hold 
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Fig. 6. The function y(x) = xsech(x)tanh(x) plotted numerically. It has 
two extremas, i.e., (xmax,ymax) = (±1.463,0.577), which are symmetric 

Fig. 5^ The limit cycle in the u - m plane, with the parameters given in about the y axis The parameter r Q is defined as the maximum va i ue f the 

Fig. [2] The dashed line represents a part of the transient. function 



Substracting Eq. ^ by Eq. ( [TO] ), we can arrive at the equation 
about the error function r](t), given by 
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dt 
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In order to calculate the evolution of the perturbation on the 
surface of section after each period T = 2tt/uj, we integrate 
over t from nT to (n + 1)T. It comes to 



\ V ((n + l)T) | =e x \r 1 (nT) \ 



where 



A 



/ S(t)dt 
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is the Floquet exponent of the linearized Poincare map. 
Let us now estimate the A without knowing the an- 
alytical solution <^(£). The even function </>^/</> m o * 
sech(0^/0 m o) tanh(0^/0 m o) in the integral is actually a 
positive bounded function in 0^, having as the minimum 
value and r as the maximum value, where r ~ 0.577 
denotes apositive real value at its two symmetric extrema 
(see Fig. ([6}); also is 0^/0 mO -sech(0^/0 mO ) tanh(<^/0 mO ) 
positive bounded. Therefore, the two parts in the integrand are 
bounded: 
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Noting that u(t) = uq cos(u;t) is alternatingly positive 
and negative in half-periods, i.e., < u(t) < uq, for t G 
[0,r/4]U[3T/4,r] and -u < u(t) < 0, for t G [T/4,3T/4], 
we integrate and estimate Eq. ( [T4| ) in three intervals with the 
relations in Eq. ( 15a) and ( 15b). This can lead to the following 
inequation: 



X < D { 



M 2 ( [ T/ \ r ) 
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Therefore, when the following condition is satisfied, 



(R^l — R^)Tq + R 




u(t)dt 



(16) 



(17) 



the Floquet exponent is guaranteed negative, which indicates 
that the system is asymptotically stable (a stable limit cycle). 
This condition is frequency independent. The left hand side 
of Eq. ([17]) is determined solely by the GMR material, while 



the right hand side is determined by the feedback coil and the 
input signal. 

IV. Discussions 

1. The spin- valve model as a standard memristor 
Let us relocate at the degenerate model (with the L- 
containing terms being non-dominant). Being both charge- 
controlled memristor models, we find it especially interesting 
that the spin- valve memristive system has a close analogy to 
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the often referred T1O2-X memristor of S. Williams, et al 0. 
To see this, we expand the sech term in Eq. ([T]) to the second 
order (valid for small fields: sech(x) = 1 — x 2 /2 + 0(x A )) and 
compare the degenerated equation set with that of the T1O2-X 
model: 



u(t) R{(j> m )i 
dln(^/^ ) = 2 , 
dt nedS 

R((/) m ) = R n - 



u(t) R(w)i 
dw _ fi v Ron 
~dt ' 



D 



(18a) 
(18b) 



{Rn ~ R n) 



2^o 



R(w) = Roff-(Roff - Ron) 



w 
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w (18c) 



Apparently, for small cj) m , taking 0^ as the state- variable, 
they have general similarities in the function relations, except 
a logarithmic function in the former. However, they have a 
complete distinction in their physical mechanisms. They differ 
also in the resistance range, being a low conductive resistance 
for the spin-valve model. Furthermore, it is notable that 
opposite to the Ti02- X memristor, the spin-valve memristor 
does not have a nonvolatile memory when switched off, as 
the current in the coil cannot hold permanently. 

2. Generalized memristive systems 

Since Di Ventra, et al, generalized the concept of memristor 
of Chua, to memristive, memcapacitive and meminductive 
systems, any newly born elementary passive memory system 
has been seeking a belonging category in the form of 



y(t) = #( x >M) • z{t) 

x = f (x, z, t) 



(19a) 
(19b) 



or of its variations, describing two complementary constitutive 
variables (current, charge, voltage or flux), where g is a 
generalized response, and f is a vector function, with y being 
the input signal, z the output signal, and x the state-variable 
vector 0. It has been successful for a very large class of 
memristive device structures |2][8]-|15]. This formal descrip- 
tion, however, is based on the linear fundamental electronic 
elements, i.e., the resistor, capacitor and inductor, which in 
common has a homogeneous characteristic that the trajectory 
must cross the origin in the u — i phase plane, due to Eq. 



(T9aj). We are curious about the necessity of holding this strong 
constraint for an arbitrary non-combinatorial passive memory 
system beyond the above. Apparently, our proposed spin- valve 
systeml has grown out of the prevailing categorization of 
memristive systems, which resembles a standard memristive 
system (pinched with self-crossing at the origin). In fact, the 
non-origin-self-crossing characteristic is permitted in a class 
of systems that are not subject to the Ohm's law. On the 
other hand, any passive system that can express a memory 
effect should be regarded as a memristive system in a general 
sense. In [26], L. Chua also concluded such property of 
resistance switching memory devices as: "If it's pinched, it's a 
memristor." This implies that the memory effect is the essential 
property of a memristor, while the origin-crossing property 
holds only for a class of memristors. In this context, it leads us 



to the generalized description of a heterogeneous memristive 
system: 

u(t) = #(x, i, t) • i(t) + <?(x, i, t) (20a) 
x = f(x,i,t) (20b) 

which does not require a solution for {i(to) = 
0,g (x(to)j 0, to) = 0}. Such systems are of special interest 
to exhibit hysteretic loop with the self-crossing knot not 
located at the origin of the u — i phase plane for the non-zero 



heterogeneous term #(x, i,t) in Eq. (20a). Heterogeneous 
memcapacitive and meminductive systems are also permitted 
for existence and can be defined in the same manner. 

This generalization can lead to a further broadened concept, 
as we have noticed that when the heterogeneous term g(x, z, t) 
gains sufficient dominancy, the dynamics of the system may 
express characteristics beyond a standard memristive system. 
This implies that intrinsic passive memory systems can possess 
more than one attributes. These attributes are, however, so 
physically entangled to each other that the boundaries of 
memristive (resistive), memcapacitive (capacitive) and memin- 
ductive (inductive) systems are vague, and hence they should 
have a compound form in general 



H(V) = 
x = P(V) 



(21a) 
(21b) 



where H and F are two vectors of implicit algebraic functions 
of a set of variables V: 



V= {u,i,u^ n \i (m) 
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Due to the shifted boundaries, we would refer to such a 
system simply as a compound memory electronic system (or 
a memtronic system in short). It is foreseeable that some com- 
plex nano-, organic or biological systems possessing electric 
memory have to be described in such a general manner. 

V. Conclusions 

We have proposed theoretically a generalized memristive 
system based on the feedback spin-valve mechanism in the 
giant magnetoresistive material with the classical Hall Effect. 
The dynamics can exhibit a pinched hysteretic loop that 
possesses the self-crossing knot not located at the origin of 
the phase plane, which is of its uniqueness with regard to 
standard memristive systems. We have also observed that 
a single-looped orbit is allowed for this system, when the 
heterogeneous terms in the governing equations gain sufficient 
dominancy. Our analysis shows that such non-origin-crossing 
dynamics is generally permitted in a class of passive memory 
systems that are not subject to Ohm's Law, which has grown 
out of the prevailing homogeneous categorization of memris- 
tive systems. We have hence introduced a broadened concept 
of the heterogeneous memristive systems and ultimately the 
compound memory electronic systems. For the stability of the 
proposed model, we have provided a sufficient condition based 
on an estimation of the Floquet exponent in the linearized 
Poincare map. 
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